Surface States of the Topological Insulator Bii_j; Sb^ 



OO 
O 

o 

(N 

< 






c^ 



Jeffrey C.Y. Teo, Liang Fu and C.L. Kane 
Dept. of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104 

We study the electronic surface states of the semiconducting alloy bismuth antimony (Bii_a;Sba;). 
Using a phenomenological tight binding model we show that the Fermi surface for the 111 surface 
states encloses an odd number of time reversal invariant momenta (TRIM) in the surface Brillouin 
zone. This confirms that the alloy is a strong topological insulator in the (1; 111) Z2 topological 
class. We go on to develop general arguments which show that spatial symmetries lead to additional 
topological structure of the bulk energy bands, and impose further constraints on the surface band 
structure. Inversion symmetric band structures are characterized by 8 Z2 "parity invariants", which 
include the 4 Z2 invariants defined by time reversal symmetry. The extra invariants determine the 
"surface fermion parity", which specifies which surface TRIM are enclosed by an odd number of 
electron or hole pockets. We provide a simple proof of this result, which provides a direct link 
between the surface state structure and the parity eigenvalues characterizing the bulk. Using this 
result we make specific predictions for the surface state structure for several faces of Bii-iSbx. 
We next show that mirror invariant band structures are characterized by an integer "mirror Chern 
number", hm, which further constrains the surface states. We show that the sign of hm in the 
topological insulator phase of Bii_2,Sba: is related to a previously unexplored Z2 parameter in the 
L point k • p theory of pure bismuth, which we refer to as the "mirror chirality", rj. The value 
of 77 predicted by the tight binding model for bismuth disagrees with the value predicted by a 
more fundamental pseudopotential calculation. This explains a subtle disagreement between our 
tight binding surface state calculation and previous first principles calculations of the surface states 
of bismuth. This suggests that the tight binding parameters in the Liu Allen model of bismuth 
need to be reconsidered. Implications for existing and future angle resolve photoemission (ARPES) 
experiments and spin polarized ARPES experiments will be discussed. 

PACS numbers: 73.20.-r, 73.43.-f, 73.61.Le 
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I. INTRODUCTION 

A topological insulator is a material with a bulk elec- 
tronic excitation gap generated by the spin orbit inter- 
action, which is topologically distinct from an ordinary 



insulator- 



j^^^^^^£ 



This distinction, characterized by 



a Z2 topological invariant, necessitates the existence of 
gapless electronic states on the sample boundary. In two 
dimensions, the topological insulator is a quantum spin 
Hall insulatorii^i^, which is a close cousin of the integer 
quantum Hall state. The edge states predicted for this 
phase have recently been observed in transport experi- 
ments on HgCdTe quantum wells^. In three dimensions 
there are four Z2 invariants characterizing a time reversal 
invariant band structure^ii. One of these distinguishes 
a strong topological insulator, which is robust in the pres- 
ence of disorder. The strong topological insulator is pre- 
dicted to have surface states whose Fermi surface encloses 
an odd number of Dirac points and is associated with a 
Berry's phase of n. This defines a topological metal sur- 
face phase, which is predicted to have novel electronic 
propertiea2ii2iiii 

In Ref. [l3 we predicted that the semiconducting alloy 
Bii_2; Sba; is a strong topological insulator using a gen- 
eral argument based on the inversion symmetry of bulk 
crystalline Bi and Sb. The surface states of Bi have been 
studied for several years. Experimentally there are sev- 
eral photoemission studies of Bi crystals and films which 
have probed the surface state9i^iiiii^ii^iiiii^ii^i2£i2ii2^. 
There are fewer studies of Bii_a: Sh^r^, but in very re- 



cent work, Hsieh et ali^ have mapped the (111) surface 
states, and verified the topological structure predicted 
for a strong topological insulator. 

First principles calculations provide a clear picture of 
the surface state structure of B i^^i^°i^^i^^i^^ , which cap- 
tures many of the experimental features, including their 
spin structure^. For the alloy, Bii_j; Sb^,, one expects 
the surface states to evolve smoothly from Bi, at least 
for small x. The alloy presents two difficulties for these 
calculations, though. First, since these calculations must 
be done on relatively thin slabs, features near the small 
band gap are inaccessible because finite size quantization 
mixes the bulk and surface states. Moreover, describing 
the alloy would require some kind of mean field treatment 
of the substitutional disorder. 

In this paper we study the surface states of Bii_a; 
Sba; first by developing a phenomenological tight binding 
model which can be solved numerically and then by devel- 
oping general arguments that exploit spatial symmetries 
and explain a number of model independent features of 
the surface states. Our phenomenological tight binding 
model is based on an interpolation of a model developed 
by Liu and Allen^^. This model has the advantage that 
it can be solved in a semi infinite geometry, which allows 
the surface state features near the small band gap to be 
calculated. Our aim is not to perform a quantitatively 
accurate calculation of the surface states, but rather to 
provide a concrete calculation in which robust, model in- 
dependent features of the surface states can be identified 
and characterized. Here we list our main conclusions: 



(1) We find that the Fermi surface of the 111 surface of 
Bii_2; Sba; consists of an electron pocket centered around 
the r point and six elhptical hole pockets centered a point 
in between F and the M point. (Here the bar refers to 
symmetry points in the 111 surface Brillouin zone) . This 
is similar to the surface states in Bi. Unhke the aUoy, 
however, Bi has bulk states at the Fermi energy: hole 
states near f and electron states near M. 

This calculation verifies the topological structure of the 
surface states predicted in Ref. [T^- In that work we 
showed that the four Z2 invariants {i/q; 1^1 1^2 1^3) character- 
izing the valence bands of pure Bi and Sb are (0; 000) and 
(1; 111) respectively. The semiconducting alloy Bii_^ 
Sba; was argued to be in the same class as Sb, which 
is a strong topological insulator. These invariants de- 
termine the number of surface bands crossing the Fermi 
energy modulo 2 between each pair of time reversal in- 
variant momenta (TRIM) in the surface Brillouin zone. 
Specifically, it predicts that for the 111 surface an odd 
number of Fermi surface lines separate the F point from 
the three equivalent M points. This is consistent with 
both our calculation and with experiments^. 

(2) We will show that for crystals with inversion sym- 
metry there is additional topological structure in the bulk 
band structure, which further constrains the surface band 
structure. At each of the 8 TRIM, F^ in the bulk Brillioun 
zone, the product of the parity eigenvalues of the occu- 
pied bands defines a parity invariant 5{Ti), which is a 
topological invariant in the space of inversion symmetric 
Hamiltonians. The four Z2 invariants, which require only 
time reversal symmetry are determined by these 8 signs, 
and determine the number of Fermi surface lines sepa- 
rating two surface TRIM. They do not, however, specify 
which of the TRIM are inside of the surface Fermi sur- 
face and which are outside. We will show that the bulk 
parity invariants 5{Ti) provide that information. 

Specifically, for each surface TRIM we will define the 
surface fermion parity as the parity of the number of 
Fermi lines that enclose that TRIM. This distinguishes 
the TRIM that are outside the Fermi surface from those 
that are inside a (single) electron or hole pocket. For a 
crystal terminated on an inversion plane we will estab- 
lish a theorem which relates the surface fermion parity to 
the bulk parity invariants. Thus, for inversion symmetric 
crystals, the 8 bulk parity invariants provide more infor- 
mation about the surface states than just the four Z2 
invariants. We will give a simple proof of this theorem in 
appendix A, which establishes a more direct connection 
between the bulk parity eigenvalues and the surface state 
structure than that presented in Ref. 12. 

For the 111 surface of Bii-ajSbj; our general theorem 
is consistent with both our surface state calculation and 
with experiment. We will also apply this result to make 
predictions about the other surfaces of Bii-j^Sb^:. In ad- 
dition, our theorem has implications for inversion sym- 
metric crystals which are ordinary insulators. In partic- 
ular, we will show that it has non trivial implications for 
the surface states of pure Bi, whose valence band is in 



the trivial (0; 000) topological class. 

(3) In addition to inversion symmetry, the crystal lat- 
tices of Bi and Sb have a mirror symmetry. We will show 
that the presence of mirror symmetry leads to a further 
topological classification of the bulk band structure in 
terms of an integer, hm, which we refer to as a mirror 
Chern number. This integer is similar to the spin Chcrn 
number, which occurs in the quantum spin Hall effect 
when spin is conserved^, and its parity is related to the 
Z2 invariants!,. The valence band of pure Bi, which has 
the (0;000) Z2 class™, has um — 0. The semiconduct- 
ing alloy is a topological insulator with Z2 class (1; 111). 
There are two possibilities for the mirror Chern number 
um — ±1, however, which correspond to topologically 
distinct phases. We will show that the sign of n_M in the 
topological insulator phase further constrains the behav- 
ior of the surface states. 

The transition between the (0; 000) and (1; 111) classes 
in Bii_j,Sb2, occurs for small x ~ .03 because pure Bi is 
very close to a band inversion transition where the Lg 
valence band and La conduction band cross. The k • p 
theory of these states has been studied extensively in the 
literature^Sii^SiHi^^ and has the form of a nearly massless 
three dimensional Dirac point. We will show that the 
change Aum in the mirror Chern number at the band 
inversion transition is determined by a previously unex- 
plored parameter in that theory: a sign ry = ±1 which 
we will refer to as the mirror chirality. rj is related to the 
sign of the g factor, which relates the magnetic moment 
to the angular momentum in a particular direction. For 
?7 = -1-1 the g factor is like that of a free electron, while 
for 77 = —1 it is anomalous. 

We will use this result to interpret our surface state 
calculation and to provide guidance for how 77 can be 
measured. In addition to the Dirac point enclosed by 
the surface Fermi surface at F, our tight binding surface 
band calculations for both pure Bi and Bii-^jShj; predict 
that the 6 hole pockets also enclose Dirac points which 
reside at points along the line between f and M. Unlike 
the Dirac points at the surface TRIM, the degeneracy 
at these Dirac points is not protected by time reversal 
symmetry, but rather by mirror symmetry. This predic- 
tion is inconsistent with first principles calculations of 
the surface states in B i^^'^^ , which do not find a band 
crossing inside the hole pocket. Since the Dirac point 
occurs above the Fermi energy ARPES experiments do 
not directly probe this issue. Nonetheless, spin resolved 
ARPES experiments on Bi provide evidence that the sur- 
face band structure is consistent with the first principles 
calculations^^. 

We will show that this inconsistency can be traced 
to the mirror chirality and the mirror Chern number. 
The mirror chirality in the topological insulator phase 
of Bii_a;Sba; can be determined from the structure of the 
k-p perturbation theory of the energy bands in the vicin- 
ity of the L point in pure Bi. We find the Liu-Allen 
model predicts that um = +1- This value implies that 
the surface state bands in the alloy cross in such a way as 



to establish the presence of the Dirac points in the hole 
pockets in agreement with our surface state calculation. 
In contrast, we find that an earlier, but more fundamen- 
tal pseudopotential calculation by Golin'^'' predicts that 
riM = —1. This value predicts that the bands do not 
cross, and that there are no extra Dirac points, which 
is consistent with the presently available experimental 
results as well as first principles calculation s^^' --. The 
Liu Allen tight binding parameters were chosen to repro- 
duce the energy of the bands computed using first prin- 
ciples calculations, incorporating available experimental 
constraints. Therefore, there is no reason to expect that 
it gets nj^ right. We conclude that the inconsistency 
in our surface state calculation is an artifact of the Liu 
Allen tight binding model, which culd be corrected with 
a suitable choice of new parameters. 

The outline of the paper is as follows. In section II 
we will review the salient features of bulk Bii_a; Sb^; and 
describe our phenomenological tight binding model. In 
section III we will describe our surface state calculations 
for Bii_2; Sba;. In section IV we will establish the rela- 
tionship between the surface fermion parity and the bulk 
parity eigenvalues and use that result to analyze the sur- 
faces of Bii_a; Sba;. In section V we will discuss the mirror 
Chern number, and show it is related to the mirror chi- 
rality of the k • p theory of pure Bi. In section VI we will 
conclude with a discussion of the relevance of our results 
to existing and future experiments. Finally, in appendix 
A we provide a simple proof of the theorem relating the 
surface fermion parity to the bulk parity eigenvalues. 



II. BULK BIi_^ SB^ 
A. Introduction 

Bismuth and Antimony are group V semimetals. They 
have the rhombohedral A7 structure shown in Fig. [Ija), 
which can be viewed as a distorted simple cubic lattice in 
which the triangular (111) lattice planes (which we will 
refer to as monolayers) are paired to form bilayers. The 
trivalent s^p^ atoms tend to form strong covalent bonds 
directed to the three nearest neighbors within a bilayer. 
Different bilayers are more weakly coupled. The primi- 
tive unit cell consists of two atoms in different monolay- 
ers, and each bilayer has a structure similar to a honey- 
comb lattice. The Brillouin zone for this lattice is shown 
in Fig. [IKb). It contains 8 special points which are in- 
variant under inversion and time reversal, denoted by F, 
T and 3 equivalent L and X points. 

Both Bi and Sb have a finite direct energy gap through- 
out the Brillouin zone, but they have a negative indirect 
gap. In Bi the conduction band minimum at L is be- 
low the valence band maximum at T, which gives rise 
to an anisotropic hole pocket and three electron pockets 
with small effective masses'^^. At L the conduction band 
minimum, which has even parity L^ symmetry, nearly 
touches the valence band maximum, with odd parity La 




FIG. 1: (a) Crystal structure of Bi. (b) 3D Brillouin zone 
and its projection onto the (111) surface. Also displayed is 
the choice of coordinate system throughout the paper: z is 
along the (111) direction, y is along the F to M direction, 
and O is a center of inversion. 



symmetry, forming a three dimensional Dirac point with 
a small mass gap Eg Ki 11 meV. In Sb, the conduction 
band minimum at L has La symmetry, and is below the 
valence band maximum at the lower symmetry H point. 
The alloy Bii_2, Sbj; retains the rhombohedral A7 
crystal structure. The evolution of its band structure 
has been studied experimentally^!^. As x is increased 
from zero two things happen. First, the small gap at 
L closes and then reopens. The Ls and La bands switch 
places, and the mass of the three dimensional Dirac point 
changes sign. Second, the top of the valence band at 
T descends below the bottom of the conduction band, 
resulting in a semimetal-semiconductor transition. For 
.09 < a; < .18 the alloy is a direct gap semiconductor a 
gap of order 30 meV at the L points. 



B. Topological Invariants 



Time reversal invariant band structures are classified 
topologically by four Z2 invariants^i^. In Ref. |3 we 
exploited inversion symmetry to show that these four in- 
variants can be determined by the parity ^m(rj) of the 
occupied bands at the 8 TRIM F^, via the quantities 



S{r^)='[[i2n{r^), 



(2.1) 



which we will refer to as parity invariants. Here the 
product includes each Kramers pair (which satisfy ^2n = 
S.2n-i) only once. For an inversion symmetric crystal, all 
8 of the parity invariants are topological invariants in the 
following sense. If the crystal Hamiltonian is smoothly 
deformed, preserving the inversion symmetry, then the 
only way any of the (5(Fi)'s can change is if the gap at 
Ti goes to zero, so that states with opposite parity can 
be exchanged between the conduction and valence band. 
If inversion symmetry is relaxed, then the 8 invariants 





5(r) S{L) 5{T) 5{X) 


(fO;J^lJ^2!^3) 


Bismuth 


-1 -1 -1 -1 


(0;000) 


Antimony 


-1 1 -1 -1 


(1;111) 


Bii_,Sb, 


-1 1 -1 -1 


(1;111) 



TABLE I: Parity invariants 5{ri) and Z2 topological invari- 
ants (i/o; 1^11^21^3) for Bismuth, Antimony, and Bii-^Sb^ deter- 
mined from the product of parity eigenvalues ^m{Fi) at each 
bulk TRIM Ti. 



lose their meaning. However, in Ref. iTa we showed that 
provided time reversal symmetry is preserved four combi- 
nations of the S{Ti) remain robust and define the four Z2 
invariants denoted by (fo; J^ii^2^3)- The most important 
invariant, j/q, distinguishes the strong topological insu- 
lator, and survives even in the presence of disorder^i^. 
(—1)'''' is given simply by the product of all 8 6{ri). 

Pure Bi and Sb have inversion symmetry. The parity 
eigenvalues for inversion about the point O in Fig. [^a) 
are tabulated in the literature^i^i^. Based on this data 
we display S{Ti) in Table |T1 along with the predicted Z2 
invariants for pure Bi, pure antimony and the alloy. The 
valence band of pure Bi is characterized by the trivial 
class (0; 000), while antimony has the (1; 111) class. The 
difference is due to the inversion of the Ls and La bands, 
which changes the sign of 6{L). The alloy inherits its 
topological class from antimony, and is a strong topolog- 
ical insulator. 



C. Pure Bi, Sb : Liu Allen model 

Liu and Allen-?- developed a third neighbor tight bind- 
ing model for the electronic structure Bi and Sb, which 
describes the atomic s and p orbitals nearest to the Fermi 
energy. The Bloch Hamiltonian H(k) = e^^^'^He^^'' has 
the form 



Hik) 



' Hn{k) i7i2(k) 

H2l(k) H22{k). 



(2.2) 



Here Hab{k) are 8 by 8 matrices describing the coupling 
between the 2 s states and 6 p states on the a and b sub- 
lattices of the crystal. The explicit form of these matrices 
is given in Tables IX and X in the appendix of Ref. l2a . 
Hii = H22 describe the coupling within the same sub- 
lattice. These terms involve the on site energies Eg and 
Ep as well an on site spin orbit coupling, A. The closest 
neighbor on the same sublattice is the third neighbor, 
which resides in the same monolayer as the origin. The 
third neighbor hopping involves four parameters V^' with 
c = ss, spa, ppa and ppTr, describing hopping between 
the s and p states. Since further neighbor hopping is not 
included in this model, 7Jii(k) and iJ22(k) describe de- 
coupled monolayers and depend only on the momentum 
q = k|| in the plane of the monolayer. 

H12 = i?2i describes the coupling between the sublat- 
tices. These involve two terms : First neighbor hopping 



terms Vc couples atoms within the same bilayer, and sec- 
ond neighbor hopping terms V^' couple atoms in neighbor- 
ing bilayers. In the following it will be useful to separate 
these two contributions by writing k = (q, kz) 



i?i2(q,fc.)=i?l2^(q)e^'= 



i/(^)(q)e-^'=^^% 



(2.3) 



where ci and C2 are the spacing between the monolayers 
within a bilayer and between different bilayers, and q and 
kz are the momenta parallel and perpendicular to the 
surface. 11^2 and iJfj can be extracted from Table X of 
Ref. |26| by noting that they are the terms which involve 
the parameters 50 — 512 and 1713 — 526 respectively. 

The 12 hopping parameters and 3 on site parameters 
make a total of 15 parameters specifying this model. 
These were chosen to reproduce the energies predicted 
by first principles calculations, as well as details of the 
band gaps and effective mass tensors which are known 
experimentally. The values of the parameters for both 
Bi and Sb are listed in table II of Ref. l2a. 



D. Tight binding model for alloy 

In order to describe the electronic structure of the al- 
loy Bii_a;Sba;, we wish to develop a "virtual crystal" ap- 
proximation which treats the substitutional disorder in 
mean field theory and results in a translationally invari- 
ant effective Hamiltonian. Since the regime of interest 
is X ~ 0.1, the effective Hamiltonian should be close to 
that of pure Bi. The effect of small x will be to modify 
the band energies, but not drastically change the wave- 
functions. The effective Hamiltonian should reproduce 
two essential features: (1) the inversion of the Lg and 
La bands (which are nearly degenerate in pure Bi), and 
(2) the descent of the valence band at T below the con- 
duction band at L, as x is increased, which leads to the 
transition between the semimetal and the semiconductor. 

The simplest approach would be to simply interpolate 
between the tight binding parameters for Bismuth and 
Antimony. For each of the 15 tight binding parameters 
ttc, we could define 



ac{x) 



xa 



Sb 



-I- (1 -x)a 



Bi 



(2.4) 



However, for this simple interpolation the inversion be- 
tween Lg and La occurs at a rather large value x ~ .4, 
which occurs after the semimetal-semiconductor transi- 
tion. We found that this could be corrected if each of the 
hopping terms (but not the other terms) are revised such 
that 



Vcix) = xV^'' + {1- x^)V^ 



(2.5) 



This approach is admittedly ad hoc, but it is sufficient for 
our purposes because it correctly accounts for the most 
important features of the band evolution. In Fig. [5] we 
plot the energies of the T'^r^, Ls and La as a function of 
x for this model. The qualitative behavior of the known 



E(eV) 




FIG. 2: Band evolution of interpolated tight binding model 
using the parameters in (|2.4l2.5p 



band evolution is reproduced, including the decent of the 
hole pocket at T and the inversion of the conduction and 
valence bands at L. This should not, however, be inter- 
preted as a quantitative description of the band evolution 

of Bii-a; Sbj;. 



III. SURFACE STATES OF BIi_^ SB^ 

In this section we describe our calculation of the 111 
surface band structure for a semi infinite lattice in the 
half plane z < described by the interpolated tight bind- 
ing model described above. We begin with a brief discus- 
sion of our method, which is based on a transfer matrix 
scheme^, and then go on to discuss the results. 



A. Transfer matrix method 

The electronic states of a semi infinite crystal can be 
represented as (f>n,a{<i) in a basis of states which are plane 
waves with momentum q in the plane of the surface, but 
are localized on the a = 1, 2 monolayer of the nth bilayer. 
Each (pna bas 8 components associated with the 8 atomic 
orbitals. The time independent Schrodinger equation, 
written in this basis may be expressed in the form 



"+^■1 Ur(q,s) y 

On+1,2 / V '^",2 



(3.1) 



where the transfer matrix is given by T(q, E) 
tii(q,£')i22(q,-B), with 



and 



-"21 [J^ - J^22) --"21 -"21 



-"12 (-t^ ~ -till) --"12 -"12 



(3.2) 



(3.3) 



Any bulk state is an eigenstate of the 16 x 16 transfer 
matrix with unimodular eigenvalues. For E within the 
energy gap, T(q, E) has exactly eight eigenvalues with 
modulus larger than 1. These correspond to states that 
decay exponentially in the —z direction. E(q) will cor- 
respond to a surface state localized at the top surface in 
Fig. O^a) near z = provided there is a linear combi- 
nation of the decaying states which vanish on the mono- 
layer n = 0, a — I just outside the surface: 0o,i = 0. The 
surface states are thus determined by forming an 8 by 8 
matrix A/(q, E) composed of the 8 components of 0o,i for 
each of the 8 decaying states. E{q) is then determined 
by solving det[A'f (q, E)] = 0. 



B. Electronic structure of (111) surface 

Fig. [3Jc) shows the energy spectrum of the (111) sur- 
face states of Bii_2:Sba; for x = .08 calculated along the 
line connecting q = F = to Af along the +y axis us- 
ing the transfer matrix method for the interpolated tight 
binding model. Fig. Elb) shows the Fermi surface. We 
find two bound surface states within the bulk energy gap. 
Along the line Qx = these states are labeled by their 
symmetry under the mirror M.{x) which takes x to —x. 
Since the mirror operation also operates on the spin de- 
gree of freedom it is important to be specific about its 
definition. We write M{x) = PC2{x), where P is inver- 
sion, and C2{x) is a 180° counterclockwise rotation about 
the positive x axis. P does not affect the spin degree of 
freedom, but the C2 rotation does. The resulting eigen- 
values of M{x) are +i and — i, which we label as Ei and 
E2. These mirror eigenvalues are correlated with the spin 
Sx- For a free spin, eigenstates with M{x) = ±i corre- 
spond to spin eigenstates with Sx = T^/2. The surface 
states are not spin eigenstates, but on the line kx = 0, 
< ky < ky{M) the expectation value of the spin satis- 
fies {S) oc i{M{x))x oc —(+)£, for Si(2)i as indicated in 
Fig. Ha). 

The Fermi surface shown in Fig. ^ih) consists of elec- 
tron and hole pockets. A single electron pocket surrounds 
r. This Fermi surface is non degenerate, and opposite 
sides of the Fermi surface are Kramers pairs with oppo- 
site spin. The electronic states pick up a Berry's phase 
of TT when they are adiabatically transported around the 
Fermi surface. This can be understood to be a conse- 
quence of the 360° rotation of the spin going around the 
Fermi surface. The Fermi surface is thus spin filtered, in 
the sense that the spin of the electron is correlated with 
it's propagation direction, roughly satisfying {S) (x q x z 
for an electron propagating in the q direction in the plane. 
In addition, there are 6 elliptical hole pockets centered 
along the 6 lines connecting f to M . These are also non 
degenerate, though unlike the electron pocket, the time 
reverse of a hole pocket is a different hole pocket. The 
crossing of the Ei and E2 bands is protected by the mir- 
ror symmetry for qx — 0. The degeneracy will be lifted 
for finite qx, so the crossing describes a two dimensional 
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FIG. 3: (a) Geometry for our surface state calculations, which 
defines our coordinate system and specifies the spin direc- 
tions of the El and S2 bands, which have mirror eigenval- 
ues +i and —i respectively, (b) Brillouin zone for the (111) 
face of Bii_3;Sba; with the electron pocket and six hole pock- 
ets predicted by our tight binding calculation, (c) Surface 
band structure along the line between f and M predicted by 
tight binding model. The shaded regions are the bulk states 
projected to the surface, (d) Schematic illustration of exper- 
imental surface band structure and Fermi surface probed by 
angle resolved photoemission spectroscopy^. The top shows 
the Fermi surface in a slice of the Brillouin zone near kx = 0, 
and the bottom shows the surface state dispersion. Compared 
with (c), there are two additional bands near M. 



Dirac point, which is enclosed by the hole pocket. 



C. Comparison with topological predictions 

A single band of surface states connects the valence and 
conduction band between f and M in Fig. Efd). This 
confirms the topological predictions for the connectivity 
of the surface state bands. In Ref. we showed that the 
number of times A7V(Aa, Ab) the surface states intersect 
the Fermi energy between two surface TRIM A^ and A^ 
satisfies 



(-1) 



AN{Aa.Ab) _ 



7r(Aa)7r(A6), 



where 



7r{Aa) = {-irS{Tal)5{Ta2). 



(3.4) 



(3.5) 



Here Tai and ra2 are the two bulk TRIM which project 
to the surface TRIM A^. The 8 parity invariants S{ri), 
defined in Eq. 12.11 are products of parity eigenvalues. 
This definition of it( 
troduced in Refs. 



irlAa) differs slightly from the one in- 
p7l[l3 because of the additional factor 



(— 1)"'. nfc is the number of occupied Kramers degener- 
ate pairs of energy bands, which is equal to the number 
of terms in the product of Eq. 12.11 For Bii-a^Sb^;, rib — 5. 
This factor does not affect AA^(Aa, Ab) in (|3.4p . However, 
this modification simplifies our further results, discussed 
below. 

For 7r(Aa)7r(Ab) = —1 there will be an odd number of 
crossings between A^ and Ab, guaranteeing the presence 
of the gapless surface states. In the appendix we will 
provide a new derivation of this connection between the 
surface states and the bulk parity eigenvalues which is 
simpler and more direct than our previous proof-^. This 
will show that with inversion symmetry the eight parity 
invariants S{Ti) contain more information about the sur- 
face state structure than just the number of crossings, a 
fact we will exploit in section IV to make general predic- 
tions about the locations of electron and hole pockets in 
the surface Brillouin zone. 

From Fig. [1] Tableland Eq. 13.51 it can be seen that 
for the alloy. 



7r(f ) = -S{T)S{T) = -1, 
7r(M) = -S{X)S{L) ^ +1. 



(3.6) 
(3.7) 



This predicts that there should be an odd number of 
crossings between F and M , which is confirmed both by 
our explicit calculation and, as we will discuss below, by 
experiment. 



D. Comparison with experiment 

Before comparing our calculation to experiment and 
other calculations, it is worthwhile to discuss what our 
calculation does not include. In addition to our approx- 
imate treatment of the alloy's bulk electronic structure, 
we have made no attempt to self consistently describe the 
potential near the surface. This will be modified by relax- 
ation of the bonds near the surface. More importantly, 
the population of the surface states determines the elec- 
tric charge distribution near the surface, which leads to 
Hartree and exchange contributions to the potential. We 
assume that the surface is electrically neutral. We will 
argue in section HI that this means that the area of the 
electron pocket is equal to the total area of the six hole 
pockets. However the potential due to a surface dipole 
layer is not included in our calculation. The effect of 
such a surface potential will be to modify the energies 
of the bands and perhaps to split off additional surface 
state bands from the continuum. However, the topolog- 
ical connectivity of the surface state bands will not be 
altered. 

In their recent ARPES experimeni*2i, Hsieh et al. mea- 
sured the spectrum of Bi.g Sb.i (111) surface states below 
Ep between F and M. The observed spectrum, which we 
have sketched schematically in Fig. [3)^d), resembles Fig. 
IHl^c), though there are some important differences. As 
in Fig. [3]^c), two surface state bands emerge from the 
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FIG. 4: (a) Bi surface states between f and M calculated 
using tight binding model, (b) Schematic picture of Bi bands 
in (a) in which Hartree effects raise the bands to accommodate 
charge neutrality. The crossing of Ei and E2 results in a 
Dirac point enclosed by a hole pocket, (c) Schematic picture 
without the crossing between Ei and E2, which resembles a 
first-principle calculation of surface states in B i^'^'^^ 



bulk valence band near f . The first intersects the Fermi 
energy forming the electron pocket centered on F, while 
the second intersects the Fermi energy forming a hole 
pocket. A third band crosses Ep from above, forming 
the opposite side of the hole pocket, and merges with the 
bulk valence band near M. Unlike our calculation, the 
observed spectrum includes an additional electron pocket 
near M. A Kramers degenerate pair of surface states is 
found in the gap at M. Away from M these states split 
to form two surface bands, which both cross Ep near the 
end of the hole pocket. Thus there are a total of five 
bands crossing Ep between F and Af , which is consistent 
with the prediction for a (1;111) topological insulator. 
The discrepancy between our calculation and the exper- 
iment is most likely a consequence of our neglect of the 
self consistent surface potential, which could lead to a 
Kramers pair of bound states to be split off from the 
conduction band at M . 

It is also instructive to compare our calculation with 
previous experimental and theoretical results for pure Bi. 
In Fig. m^a) we show the surface state spectrum for pure 
Bi calculated using the transfer matrix method for the 
Liu-Allen tight binding model. The number of band 
crossings is consistent with the trivial (0; 000) topolog- 
ical structure of the Bi valence band. Since the Fermi 
energy of semimetallic Bi is fixed by the bulk, our calcu- 
lated surface states violate surface charge neutrality: too 
many surface states are occupied, so the surface will have 
a negative charge. Hartree effects will push the surface 
states up in energy, but they will not alter the topological 
connectivity of the surface states. This allows us deduce 
qualitative conclusions from the calculation. 

First, as in our alloy calculation, two surface bands 
emerge from the bulk valence band near f . These are also 
seen in photoemission experiments as well as first princi- 
ples calculations on pure Bi^-. Moreover, the spin (S'x) 
of those surface states has been both calculated and mea- 
sured using spin polarized ARPESSii. We have checked 



that the spin direction predicted by our tight binding 
calculation for each of these bands agrees with the ex- 
perimental and first principles theory results. Thus, the 
behaviour near F, including the ordering in which the E2 
emerges first and forms the electron pocket appears to 
be robust, with all calculations in agreement with each 
other and with experiment. 

There is a discrepancy, however, between the 
tight binding calculation and the first principles 
calculation22i24. The crossing between the Ei and S2 
bands in Figs. [3l^c) and HJa) is not found in the first 
principles calculation. Since it is likely that this crossing 
would be pushed above the Fermi energy by Hartree cor- 
rections (so that the crossing occurs inside a hole pocket), 
the tight binding model predicts that the hole pockets of 
Bi (111) enclose a Dirac point, as shown schematically 
in Fig. mb). The existence of this band crossing is not 
directly probed by ARPES which only probes occupied 
states, though it could be probed using inverse photoe- 
mission. There is, however, indirect experimental evi- 
dence that the crossing does not occur. Spin polarized 
ARPES measurements^^, have measured the spin on both 
sides of the hole pocket. Though the signal appears weak, 
the sign of the spin is resolved, and determined to be the 
same on both sides, indicating that there is no crossing, 
as shown schematically in Fig. HJc). This agrees with the 
predictions of the first principles calculations that both 
sides are in the same Si band. In contrast, our tight 
binding model predicts that the opposite sides of the hole 
pockets correspond to the Si and E2 bands, which have 
opposite spin. 

It thus appears likely that the prediction of the level 
crossing which implies that the hole pockets enclose a 
Dirac point is an artifact of the tight binding model. 
This brings into question the related prediction of the 
tight binding model that the hole pockets of the alloy 
also enclose a Dirac point. In Section V we will argue 
that this artifact is a consequence of a subtle error in the 
Liu- Allen tight binding model. 



IV. INVERSION SYMMETRY AND THE 
SURFACE FERMION PARITY 

An inversion symmetric crystal can have no bulk elec- 
tric polarization. In this section we show that this fact in 
combination with surface charge neutrality has non triv- 
ial implications for the surface state structure because 
it allows the outside of the surface Fermi surface to be 
unambiguously defined. It is then possible to define elec- 
tron pockets to be regions in the surface Brillouin zone 
where an extra band is occupied and hole pockets as re- 
gions where an otherwise occupied band is empty. Charge 
neutrality dictates that the area of the electron pockets 
should equal that of the hole pockets. We will show that 
the locations of the electron and hole pockets in the sur- 
face Brillouin zone are topologically constrained by the 
bulk parity invariants (S(Fi). In addition to fixing the 



number of Fermi energy crossings, we find that 5{Ti) de- 
termine which TRIM are on the inside of an electron or 
hole pocket and which TRIM are on the outside. We de- 
fine the surface fermion parity, which specifies whether a 
given surface TRIM is enclosed by an even or odd number 
of Fermi lines. We will begin with a general discussion 
of the relationship between the surface fermion parity to 
the bulk parity invariants. We will then apply our general 
result to the surfaces of Bii_j:Sbj; and Bi. 



A. Surface fermion Parity 

The total surface charge density may be expressed as 
a sum over the surface Brillouin zone (SBZ), 



SBZ 






Niq), 



(4.1) 



where the surface fermion number A^(q) represents the 
excess charge in the vicinity of the surface due to states 
with momentum q in the plane of the surface. If we 
assume that the bulk Fermi energy is inside the gap, 
then there will be two contributions, A^(q) = A^buikCq) + 
^surface (q)- A'^3uifaco(q) IS an integer which counts the 
occupied discrete surface states inside the energy gap. 
iVbuik(q) is the total surface charge in the continuum va- 
lence band states. For a crystal with inversion symmetry 
there can be no bulk electric polarization, and Afbuik(q) 
will also be quantized. In Appendix A we will show that 
it must be an integer—. 

The integer values of A^(q) allow us to unambiguously 
define the "outside" of the surface Fermi surface to be the 
region for which A^(q) = 0. A^(q) = +(— )1 define elec- 
tron (hole) pockets. A^(q) — +(— )2 is a double electron 
(hole) pocket, and so on. From (|4.ip . charge neutrality 
implies that the total area of the electron pockets equals 
that of the hole pocket, provide the double pockets are 
appropriately counted. 

Kramers' theorem requires that the surface states be 
two fold degenerate at the TRIM q = A^ in the surface 
Brillouin zone. Provided the Fermi energy is not exactly 
at the degeneracy point this means that A"sui.face(Ao) is 
even, so that the parity of N{Aa) is equal to the par- 
ity of A'buik(Aa). In Appendix A we will show that the 
surface fermion parity is determined by the bulk parity 
invariants, 

(_i)A'(A„) = ^(^^) ^ (-i)"'<5(r,i)<5(r<,2). (4.2) 



Eq. 14.21 determines whether the TRIM Aa is enclosed 
by a single (or odd number) of Fermi lines, or whether 
it is outside the Fermi surface (or enclosed by an even 
number). In the special case that the Fermi energy is 
exactly at a Dirac point at Aa, Aa should be interpreted 
to be inside an electron (or hole) pocket with vanishing 
size. 

Eq. 14.21 is a new result which provides information 
about the structure of the surface Fermi surface beyond 
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FIG. 5: Two inequivalent inversion centers c and c' in an 
inversion-symmetric crystal, which differ by half a lattice vec- 
tor. The parity eigenvalues of Bloch state at momentum 
k — tt/R with inversion center chosen at c and c' are different. 
Crystals terminated at c and c' will have surface charges that 
differ by an odd integer. 



that determined by the Z2 invariants (i^o; J^ii'2^'3)- We 
will show below that this result can have non trivial con- 
sequences even in materials which are not topological in- 
sulators. For example, we will see that (|4.2p constrains 
the surface states of pure Bi. 

In order to apply (|4.2p it is essential to use the par- 
ity eigenvalues associated with an inversion center in the 
plane on which the crystal is terminated. As a simple 
example Fig. [5] shows a one dimensional inversion sym- 
metric lattice, which has two distinct inversion points. In 
general, a three dimensional inversion symmetric crystal 
has 8 distinct inversion centers, which are related to each 
other by half a, Bravais lattice vector: c' = c -I- R/2. The 
parity eigenvalues associated with inversion center c' will 
be related to those associated with c by 



C(r.)-^m(r.)e' 



■Ti-R 



±^m(r.). 



(4.3) 



An inversion plane will contain four of those points. 
For a given surface orientation there are two distinct par- 
allel inversion planes. For a surface terminated on one of 
those inversion planes, 7r(Aa) does not depend on which 
of the four inversion centers within the inversion plane 
are used. This can be seen by noting that 



7r'(Aa) = 7r(Aa) exp[mf,(rai - ra2) • R], 



(4.4) 



where rib is the number of occupied bands. When c and 
c' are in the plane of the surface the dot product in the 
exponent is zero. Crystals terminated on inequivalent 
inversion planes, however will have different N(Aa). For 
odd ni,, 7T'{Aa) = — 7r(Aa), so that the parity of A^(Aq) 
changes at all four Aq. Thus, changing the inversion 
plane amounts to filling (or emptying) a single surface 
band throughout the surface Brillouin zone. Since N{Aa) 
depends on how the crystal is terminated, it is not a bulk 
property. However, AN{Aa,Ab) = N{Aa)-N{Ab) mod 2 
is a bulk property, which is determined by the Z2 invari- 
ants {vq; 1^11^21^3) . 



B. Application to Bii_j,Sba; 

We now apply our general result to Bii_3;Sba; surfaces. 
In order to apply (|4.2p it is necessary to identify the ap- 
propriate inversion centers. The 8 inversion centers of the 



Face 


Cj 


Aa = (raira2) 


7rBiSb(Aa) 


7rBi(Aa) 


(111) 


Co Cl2 


f = (FT) 


-1 


-1 




Cl3 C23 


3M = (LX) 


+ 1 


-1 


(111)' 


Cl C2 


f = (rr) 


+ 1 


+ 1 




Ca Ci23 


3M = (LX) 


-1 


+ 1 


(110) 


Co C3 


f = (rx) 


-1 


-1 




Cl2 Cl23 


Xi = (ii) 


-1 


-1 






X2 = (LT) 


+ 1 


-1 






M = (XX) 


-1 


-1 


(100) 


Cl Cl3 


f = (FL) 


-1 


+ 1 




C23 Ci23 


M = (TX) 


+ 1 


+ 1 






2M' = (LX) 


-1 


+ 1 



TABLE IL For each crystal face (hkl) we list the 4 inver- 
sion centers Cj on the cleavage plane along with the projec- 
tions relating the 4 surface TRIM Aa to the bulk TRIM Fai,2. 
For each A we list the surface fermion parity, 7r(Aa) for both 
Bii-iiSbi, and Bi. 7r(Aa) is a product of parity invariants at 

Fal,2- 



rhombohedral A7 lattice are (1) Cq — 0, the origin in Fig. 
1, which is between two bilayers. (2-4) Cj=i_2,3 = slj/2. 
Here a,,- are the three rhombohedral primitive Bravais lat- 
tice vectors, which connect an atom to the nearest three 
atoms on the same sublattice of the neighboring bilayer— . 
These points are at the center of a nearest neighbor bond 
in the middle of a bilayer. (5-7) Cy = {a.i+aj)/2 for i ^ j- 
These three points are at the center of a 2nd neighbor 
bond between two bilayers. (8) C123 = (ai + a2 + a3)/2, 
which is directly above the origin in Fig. [Jl in the middle 
of a bilayer. For a given surface orientation, these inver- 
sion centers are divided into two groups of four, which 
reside in two possible cleavage planes. 

In Ref. [H, the (111), (110) and (100) faces of Bi are 
discussed, where the Miller indices (mno) refer to the 
rhombohedral reciprocal lattice vector mhi + nh2 + oh^ 
with a.; • hj = 2TrSij . In these cases the preferred cleavage 
plane is the one which minimizes the number of broken 
first neighbor bonds. In table|TT]we list the four inversion 
centers in the cleavage plane for each of these faces. For 
comparison, we have also included the (111)' face, which 
is terminated in the middle of a bilayer (breaking three 
nearest neighbor bonds) . Table [Til also shows how the 
bulk TRIM project onto the surface TRIM, using the 
notation A^ = (rQira2)- This data, combined with Table 
im is sufficient to determine the surface fermion parity 
7r(Aa) for both the alloy Bii_a;Sb2: (BiSb) and pure Bi 
for each surface as shown in Table [ill 

First consider the 111 surface. The parity eigenvalues 
quoted in the literature, which determined Eq. 12.11 in 
Table [H are with respect to an inversion center between 
two bilayers (point O in Fig. [ija)). Thus, for a crystal 
cleaved between two bilayers, N{Aa) can be deduced by 
combining Eq. 14.21 with 
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FIG. 6: Schematic diagram showing which surface TRIM are 
enclosed by an odd number of electron or hole pockets for dif- 
ferent faces of Bii_a;Sbi predicted by the surface fermion par- 
ity in TableHI] (a-d) show the (111), (111)', (110) and (100) 
faces. The (111)' surface is a hypothetical surface cleaved in 
the middle of a bilayer. 



as shown in Table [TTl This implies the surface Fermi 
surface encloses F, but not M, as shown schematically 
in Fig. [HJa). Eq. 14.21 savs nothing about either the hole 
pockets seen in experiment and our calculation or the 
double electron pocket at M observed in experiment- 
on Bii_2;Sba; but not our calculation. In order for the 
surface to be neutral, however, the Fermi energy must 
either be at a Dirac point at f (so that the Fermi surface 
has vanishing area) or there must also be compensating 
electron/hole pockets elsewhere in the surface Brillouin 
zone (but not enclosing M). 

It is also instructive to first consider a (111)' face 
cleaved between the monolayers in a bilayer, despite the 
fact that such a surface would likely be unstable. Since 
the origin Cq is not in the cleavage plane, the parity eigen- 
values in (|2.ip need to be modified using 14.41 This has 
the effect of changing the sign of all of the 7r(Aa), so that 



^'(A,) = +S{Tal)6{ra2). 



(4.6) 



7r(A„) = ~SiTal)S{Ta2) 



(4.5) 



From Table [TTl we thus conclude that the three M points 
are enclosed by the Fermi surface, but not F, as shown 
in Fig. M.'b). 

For the 110 surface the cleavage plane with one broken 
bond includes the origin Cq. Thus 7r(AQ) can be deter- 
mined with (|4.5p along with the projections of the bulk 
TRIM shown in Table [TTl This leads to the predictions 
for the surface Fermi surface shown in Fig. ^c). Ex- 
perimental data for this face of Bii_a;Sba; is currently 
unavailable. However, it is instructive to compare this 
prediction with experiments on pure Bi. In Ref. [l5l . sin- 
gle hole pockets are clearly seen at F and M, and at Xi 
single surface Dirac point is present inside the bulk gap. 
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The situation at X2 is obscured due to the overlap of the 
bulk conduction and valence bands at L and T. 

For the 100 surface the cleavage plane with one broken 
bond does not include Cq. Thus, as was the case for the 
(111)' surface, the surface fermion parity follows from 
(|4.6p . The surface Brillouin zone shown in Fig. [S]^d) 
has TRIM F, M and two equivalent M' . Again, there is 
presently no data for this surface of Bii xSbr. The (100) 
face of pure Bi is discussed in Ref. [20, and appears to 
be consistent with the prediction of Table ITTl that none of 
the TRIM are enclosed by a Fermi surface. 



V. MIRROR CHERN NUMBER AND THE 

MIRROR CHIRALITY AT THE L POINT OF 

BISMUTH 

In this section we will explore the consequences of mir- 
ror symmetry on the band structure of Bi and Bii-^jSh^:. 
This will address the disagreement between our calcula- 
tion of the surface band structure and previous experi- 
mental and theoretical results. As discussed in Section 
3, the tight binding model predicts that the hole pockets 
enclose Dirac points, while experiment and first princi- 
ples calculations suggest that they do not. Here we will 
show that the presence of this crossing probes a funda- 
mental, but previously unexplored, property of the bulk 
electronic structure of Bi. 

We will begin by pointing out that the mirror symme- 
try of the rhombohedral A7 structure leads to an addi- 
tional topological structure of the energy bands which we 
refer to as a mirror Chern number. We will then show 
that the value of this integer in the topological insulator 
phase depends on the structure of the nearly degenerate 
Ls and La bands in pure Bi. We will identify a previously 
unexplored parameter in the k • p theory of Bi, which we 
refer to as the mirror chirality. We will show that the 
mirror chirality at the L point in Bi determines the value 
of the mirror Chern number in the topological insulator 
phase of Bii-a^Sba,. 

We find that the value of the mirror chirality predicted 
by the Liu Allen tight binding model^^ disagrees with the 
value predicted by a more fundamental calculation by 
Golin"^^. This, combined with the disagreement with the 
surface state experiments and first principles calculations 
suggests that the Liu Allen tight binding model has a 
subtle, but topological, error. 



A. The Mirror Chern Number 

The Dirac points in the hole pockets in our tight bind- 
ing calculation arise because the Si and S2 bands cross 
on the line connecting F and M in Fig. ^h). This 
crossing is protected by the invariance of the Hamilto- 
nian under the mirror operation M{x) = PC2{x) which 
takes X to —x. l]i(I]2) transform under different repre- 
sentations oi M.{x) with eigenvalues -\-i{—i). This mirror 



symmetry implies that all the bulk electronic states in the 
plane k^ — Q can be labeled with a mirror eigenvalue ±«. 
Within this two dimensional plane in momentum space, 
the occupied energy bands for each mirror eigenvalue will 
be associated with a Chern invariant n±i . Time reversal 
symmetry requires that n^i + n^i = 0, but the difference 
defines a non trivial mirror Chern number 



njvi = {n+i - n-0/2. 



(5.1) 



The situation is analogous to the quantum spin Hall state 
in grapheneii^, where the conservation of spin Sz leads 
to the definition of a spin Chern numbcr'^^>, whose parity 
is related to the Z2 topological invariant. 

The mirror Chern number determines how the surface 
states connect the valence and conduction bands along 
the line g^ = between F to M . To see this, consider 
the Ai = zLi sectors independently. The bulk states with 
kx = are then analogous to a two dimensional integer 
quantum Hall state with Hall conductivity n±ie'^/h. The 
sign of n±i determines the direction of propagation of the 
edge states, which connect the valence and conduction 
bands. Thus, the sign of n_\4 determines whether the Si 
band or the S2 band connects the valence and conduction 
band between F and M (which we take to be in the +y 
direction). For um ~ +1(— 1) we find that the Si(S2) 
band crosses. 

The predictions of the tight binding model are more 
likely to be robust near F than near M, because near F 
they are not sensitive to the detailed treatment of the 
small bulk energy gap at the L point. This is supported 
by the fact that the ordering of the Si and S2 bands 
near F predicted by the tight binding model (in which 
S2 emerges first) agrees with other calculations and ex- 
periment. Given this ordering near F, the mirror Chern 
number determines whether or not the Si and S2 bands 
have to cross. Referring to Fig. ^c), if the mirror Chern 
number were to have the opposite sign, then the S2 band 
would connect to the conduction band rather than the Si 
band, and the bands would not have to cross. Pure Bi 
is very close to the transition between the (0; 000) and 
(1; 111) phases. Therefore, it is likely that the presence 
of the crossing between Si and S2 will be unaffected by 
the transition. Therefore, the sign of the nj^ in the topo- 
logical insulator phase of Bii_a;Sba; should be correlated 
with the alternatives shown in Figs. [H with um — +(~)1 
corresponding to Fig. |1|d ([3J:). 

Since the valence band of pure Bi is in the trivial 
(0; 000) topological class, pure Bi does not have surface 
states which connect the valence and conduction bands. 
Thus the mirror Chern number for the k^ — plane of 
the valence band of pure Bi is hm — 0. The transition 
to the strong topological insulator in Bii-ajSh^, occurs 
for small x because the Ls and La bands in pure Bi are 
nearly degenerate. At the transition to the topological 
insulator the two bands cross and form a three dimen- 
sional Dirac point at L. At this transition both the Z2 
topological invariants {vq; i^ii^2^3) and the mirror Chern 
number um change. The change Aum across this tran- 
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sition is an intrinsic property of this Dirac point. Thus 
the value of njn in the topological insulator phase can 
be determined by studying the properties of this Dirac 
point. Since pure Bi is very close to this transition, this 
information can be extracted from the structure of the 
k • p Hamiltonian for pure Bi in the vicinity of the L 
point. 

In the next section we will analyze the k ■ p theory and 
show that the value of Ahm predicted by the Liu Allen 
tight binding model disagrees with the value predicted 
by an earlier pseudopotcntial calculation by Golin. This 
provides evidence that the crossing of the Si and S2 
bands is an artifact of the incorrect sign of hm predicted 
by the tight binding model. 



B. k ■ p theory and the mirror chirality 

The k • p analysis of Bi near the L point has a long 
history. Originally developed by Cohen and Blount^S 
in 1960, the theory was given an particularly elegant 
formulation by Wolf|22-, who emphasized the similarity 
with the relativistic Dirac equation. This theory, and its 
refinements^i^i^ played an important role in the early 
development of band theory, and formed the framework 
for interpreting a large body of magnetic, transport and 
optical data. In this section we point out a previously un- 
explored sign which characterizes this theory : the mirror 
chirality. We show that it is this sign which determines 
the sign of hm in the topological insulator phase. 

The four relevant states at the L point are denoted 
[Lg, La] = [(^6, -Z^s), {Li, Ls)^. The two states compris- 
ing Lg and La are degenerate due to time reversal sym- 
metry. These states are distinguished by their symme- 
try under parity P (with eigenvalues [(1, 1), (— 1, — 1)]), 
under the twofold rotation C2{x) (with eigenvalues 
[{—i,i),{i,—i)]) and under the mirror M{x) — PC2{x) 
(with eigenvalues [(— i,i), (— z,z)]). We have chosen the 
unconventional order of the states to simplify the mirror 
operator, which makes the connection with the mirror 
Chern number in section VC the most transparent. In 
this basis the inversion, rotation and mirror operators 
have the direct product form. 



p = 


Tz ® 1, 


C2{X) = 


-iTz ® y.z 


M{x) ^ 


-it® Hz: 



while the time reversal operator can be chosen as 

e = il (g) ^iyK, 



(5.2) 



(5.3) 



where K is complex conjugation. pL and f are Pauli matri- 
ces operating within and between the Lg and La blocks, 
and 1 is the identity matrix. In the following we will 
simplify the notation by omitting the ® and the 1. 
To first order in k the k • p Hamiltonian has the form 



Where Eq — "im is the energy gap (positive for Bi) and 
Ha are 4x4 matrices. Invariance of 7?(k) under P and 
Q requires {Ha, P} = {Ha, 8} = 0, and invariance un- 
der M{x) requires {U^,M{x)} = [Uy^z,M] = 0. The 
allowed terms are thus 



(5.5) 



^x = tlTxfix + hTxfJ-y 
Uy = UllTxfiz + Ul2Ty 
n^ = U2lT.j,^lz + U22Ty 



where ti and u^ are real numbers. Eqs. 15.41 and 15.61 
are equivalent to the k • p theory introduced by Cohen 
and Blount^^, who expressed the Hamiltonian in terms 
of complex vectors t and u. These are related to our 
parameters via t = (<i + it2)x and u = (— un + iui2)y + 
{—U2i+iu22)z . In the following it will be useful to express 
these in terms of three complex numbers t — x ■ i and 
u^ = (y ± iz) ■ u. 

Eq. 15.41 has a simpler form when expressed in terms 
of the principle axes in both momentum space and 
spin space. We thus perform a rotation [ky ~\- ikz) = 
e'"(fc' -j- ik'^) along with a unitary transformation \'ip) = 
exp[i^2(/3 -I- 7X2)]!?/;'). These transformations have the 
effect of changing the phases t -^ te*^ and u^ -^ 
y±g-i(7±Q;)^ For appropriately chosen a, (3 and 7, t and 
=pu can be made real and positive. The Hamiltonian 
then takes the diagonal form 



H = niTz + Vik^T^fi^ + ■r]V2kyT^fiz + v^k'/ 



where 



Vl = 


= \t\ 


'7^2 = 


= (|^.+ |-K|)/2 


V3 = 


= (|^.+ | + K|)/2 



Hik) 



kylly 



kzUz, 



(5.4) 



(5.6) 



(5.7) 



Here we have defined V2 to be positive and introduced a 
previously unexplored quantity 77 = ±1, which is simply 
given by 77 = sgn(det[uy]). 77 is a mirror chirality, which 
distinguishes two topologically distinct classes of Dirac 
Hamiltonians. 

For a system with full rotational symmetry, 77 must 
be equal to +1. This can be seen by noting that the 
twofold rotation operator specifies the generator of con- 
tinuous rotations about x via C2{x) = exp[— zttS'j:]. Since 
C2ix) = —i^zTz, this implies Sx = [IzTzI^- When 
77 == — 1 , Eq. 15.61 is not invariant under continuous rota- 
tions generates by Sx even when V2 — V'i, since the spin 
and orbital degrees of freedom are rotated in opposite 
directions. The twofold rotational symmetry, however, 
remains intact. 77 = +1 corresponds to the behavior of a 
free electron and should be considered normal behavior. 
77 = —1 is anomalous. 

The sign of 77 is not ordinarily discussed in the k • p 
theory of Bi because it has no effect on the electronic 
dispersion £^(k), which depends only on \va\- t\ does, 
however, have a subtle effect in the presence of a mag- 
netic field. A magnetic field in the x direction leads to a 
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splitting of states according to their spin angular momen- 
tum Sx, which can be defined as above in terms of the 
twofold rotation operator C2{x). This defines a magnetic 
moment, which symmetry restricts to be either parallel 
or antiparallel to x. The form of this magnetic moment is 
discussed in Refs. I30ll3ll . and it is straightforward to show 
that fl ex rjSxX. This means the 77 determines the sign 
of the g factor, which describes the relation between the 
magnetic moment and angular momentum. For r/ = +1 
the sign is the same as that for a free electron, while 
•q — —1 the sign is opposite. 

Unfortunately, this sign is difficult to probe experimen- 
tally. In addition to complications which arise due to the 
presence of three equivalent L points, measurement of the 
sign requires measurement of the spin angular momen- 
tum in addition to the change in energy with magnetic 
field. The selection rules discussed in Ref. [30 are unaf- 
fected by the sign. We are not aware of any experiments 
on Bi which directly probe this sign. 



C. Relation betAveen mirror chirality, mirror Chern 
number and surface states 

We will now argue that the sign of 77 determines the 
sign of the mirror Chern number in the topological insu- 
lator phase of Bii_3,Sba,. This leads to an experimentally 
testable prediction regarding the crossing of the surface 
states. Thus, probing the surface states of the topologi- 
cal insulator may well be the best experimental method 
for determining this fundamental parameter of the k • p 
theory of Bi. 

The connection between -q and the mirror Chern num- 
ber can be established by considering the mirror plane 
fcj, = 0. H then decouples into two independent two band 
Hamiltonians for Ad{x) = —ifiz = i* with the form 



h — rriTz + sv2k' Tx + vsk'^Ty. 



(5.8) 



where s = rjHz- m 

the Chern number n^ 

from negative to positive 

change in the mirror Chern number. 



= describes a transition where 
i^^ changes. When m changes sign 



A 



n^ 



Iflz 



v^^z■ 



Thus, the 



Aum = nM{m > 0) — nj^iin < 0) = —r], 



(5.9) 



depends on the mirror chirality 77. Since jim — for Bi 
(with 771 > 0), we conclude that the topological insulator, 
with m < has 



The dispersion for the surface states on the top surface 
along ga; = is 



E{qy) = -VfJ'zmy 



(5.11) 



with 77 > 0. This means that the Ei band, which has 
fiz — ^1, propagates in the +riy direction, while the S2 
band, with fiz — +1 propagates in the —r]y direction. 
Therefore, the surface state connecting the valence band 
to the conduction band which has the positive velocity in 
the y direction will be Si for 77 = +1 and S2 for 77 = —1. 



D. Comparison of tight binding and 
pseudopotential models with experiment 

In this subsection we show that the value of 77 pre- 
dicted by the Liu Allen tight binding model^^ disagrees 
with that predicted by an early calculation by Golin^. 
Specifically, we find that the Liu Allen model predicts 
the conventional value, 77 = 1, while the Golin model pre- 
dicts the anomalous value 77 = —1. We will then argue 
that the value of 77 can be extracted from the structure 
of the surface state spectrum. The presently available 
spin polarized ARPES data on the Bi 111 surface^^ pro- 
vides indirect evidence that the mirror chirality has the 
anomalous value 77 = —1. 

The k • p parameters can be determined by evaluating 
the matrix elements 



njf' - {L,\Va\Lj)\k^L. 



(5.12) 



where v — Vk'W(k)|k=L is determined by the Bloch 
Hamiltonian Ti(k). From this it follows that 



t 

,± 



n 



57 



K' T *n 



67 



(5.13) 

(5.14) 



These matrix elements are listed in table II of Golin's 
paper'^^ (the relevant band is j — j' = 3). They may also 
be extracted from the Liu Allen tight binding model. In 
Table Hill we compare the values of ui, V2, 773 and 77 com- 
puted from these matrix elements. The signs of 77 pre- 
dicted by the two theories disagree. Since the parameters 
of the Liu Allen model were simply fit to reproduce the 
energies of the bands, there is no reason to expect that 
it gets 77 right. In contrast, Golin's calculation, which is 
based on a pseudopotential approach, starts from more 
fundamental premise. 



UM = fh 



(5.10) 



um determines the direction of propagation of the Si 
and S2 surface states along the mirror line q^ — 0. The 
direction of propagation of the surface states on the top 
surface which connect the valence and conduction bands 
can determined by solving (|5.8p with a z dependent mass 
777,(2) = 77isign(2:) with 777 > 0. The bound state at the 
surface has wavefunction proportional to exp(— |77t,z|/7;3). 





wi (eVA) 


V2 


V3 


V 


Golin Pseudopotential 


4.16 


1.37 


7.01 


-1 


Liu Allen Tight Binding 


5.89 


0.92 


9.67 


+1 



TABLE III: Parameters of the k ■ p theory, Eq. 15.61 ex- 
tracted from the pseudopotential model'- and the tight bind- 
ing model^ 
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In the previous section we showed that provided there 
is only a single transition between pure Bi and the topo- 
logical insulator phase of Bii-ajSh^^, the mirror chirality 
deduced from the pure Bi band structure determines the 
mirror Chcrn number in the topological insulator. This, 
in turn, determines the direction of propagation of the 
Si and S2 states along the line qx — 0. The surface 
state structure predicted by the tight binding model was 
shown in Fig. [3l^c). The crossing of the Ei band is con- 
sistent with ry = -1-1. This crossing guarantees that there 
is a Dirac point enclosed by the hole pocket. This can 
be probed either by inverse photoemission or by spin po- 
larized photoemission. In the latter case, the presence 
of the Dirac point would lead to a change in the sign 
of the spin on either side of the hole pocket. It will be 
interesting to experimentally determine this property for 
Bii_2;Sb2: using spin polarized ARPES. 

Currently available spin polarized photoemission data 
on the 111 surface of pure Bi^^ provide an indirect probe 
of 77. Hole pockets are observed along the line from f 
to M in both Bii_a;Sba; and pure Bi. Provided we make 
the plausible assumption that no additional level cross- 
ings occur near the transition to the topological insulator, 
then the presence or absence of Dirac points in the hole 
pockets should be the same on both sides of the transi- 
tion. In Ref. 22, the spin in either side of the hole pocket 
was found to point in the same direction, which indicates 
that in pure Bi, the hole pockets do not enclose a Dirac 
point. This conclusion was supported by first principles 
surface state calculations, which also find no crossing^^. 
This suggests that in the alloy, it should be the S2 band 
which connects the conduction and valence bands, which 
is consistent with ry = — 1. 

It thus appears likely that the mirror chirality in Bi 
has the anomalous sign, rj — —1. This conclusion con- 
tradicts the prediction of the tight binding model, but it 
is supported by (1) the pseudopotential band structure 
of pure Bi and (2) the observed and calculated surface 
state structure of pure Bi. Spin polarized ARPES exper- 
iments on the topological insulator Bii-^jSb^: could more 
directly determine this sign by probing the mirror Chern 
number jim- 



VI. CONCLUSION 

In this paper we have analyzed the surface state struc- 
ture of the topological insulator Bii_a;Sba;. Using a sim- 
ple tight binding model based on Liu and Allen's tight 
binding parameterization we confirmed that the surface 
states have the signature of the strong topological insu- 
lator by showing that the surface Fermi surface encloses 
an odd number of Dirac points. The tight binding model 
also predicts that the surface is semi metallic, with an 
electron pocket centered on f along with 6 hole pockets. 

Using general arguments based on inversion symme- 
try, we showed that the location of electron and hole 
pockets in the surface Brillouin zone is constrained by a 



quantity which we defined as the surface fermion parity. 
This quantity is determined by the parity invariants of 
the bulk band structure, and for a given surface it deter- 
mines which surface TRIM are enclosed by an odd num- 
ber of Fermi surface lines. This argument establishes a 
simple and direct connection between the bulk electronic 
structure and the surface electronic structure for crystals 
with inversion symmetry. Using this general principle, 
we predicted the structure of the surface states for sev- 
eral different faces of Bii_a;Sba;. For the 111 face, these 
predictions agree both with our surface state calculations 
and with experiment. It will be interesting to test these 
predictions experimentally on other faces of Bii_2;Sba;. 

Finally, we showed that the mirror symmetry present 
in the rhombohedral A7 lattice leads to additional topo- 
logical structure in the bulk energy bands. We defined an 
integer mirror Chern number nM j whose value is nonzero 
in the topological insulator phase. The sign of um deter- 
mines the direction of propagation of each of the surface 
states along the mirror plane, and thus determines which 
surface states connect the conduction and valence bands. 
We find that the crossing of the Ei band predicted by the 
tight binding model, which leads to a Dirac point in the 
hole pockets, disagrees with the natural extrapolation of 
experiments and first principles calculations on pure Bi, 
which find no Dirac point in the hole pockets. 

We traced this discrepancy to a previously unexplored 
property of the k ■ p band structure of pure Bi, which 
we defined as the mirror chirality, 77. We showed that rj 
in pure Bi determines um in the topological insulator. 
Moreover, we showed that the Liu Allen model predicts 
the conventional value 77 = -1-1, while an earlier pseu- 
dopotential calculation by Colin predicts the anomalous 
value 77 = —1. The latter value is consistent with the 
available experimental data on Bi, though the connec- 
tion is rather indirect. A more direct test would be to 
directly measure the mirror chirality um in the topo- 
logical insulator by probing the surface states with spin 
polarized ARPES. 

It would be interesting to check that the value of 77 
predicted by more accurate first principles calculations 
of Bi agrees with the pseudopotential prediction. Since 
the tight binding model was designed only to get the en- 
ergies of the bands right, there is no reason to expect 
that it would get rj right. It should be possible to come 
up with a new parameterization of the Liu Allen model 
which would have rj = —1. We expect that the sur- 
face states computed within this model would have band 
crossings which agree with experiment and first principles 
calculations, though of course a quantitative description 
of the surface states requires an accurate description of 
the surface potential. 

An important lesson to be learned from this paper is 
that in addition to time reversal symmetry, spatial sym- 
metries can play an important role in topologically con- 
straining bulk and surface band structures. Our analysis 
of these symmetries has not been exhaustive. A complete 
theory of topological hand theory, which accounts for the 
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full point group symmetry of a crystal is called for. 
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APPENDIX A: SURFACE FERMION PARITY 
FROM BULK PARITY INVARIANTS 

In this appendix we show that for an inversion and 
time reversal invariant crystal the surface fermion num- 
ber A^(q = Aa) discussed in section IV is an integer, 
whose parity is determined by the product of bulk parity 
invariants S{Tai,2), which are products of parity eigenval- 
ues given in Eqs. I2.1l and l3.5l The simple proof outlined 
here provides a direct connection between the topological 
structure of the surface states and the parity eigenvalues 
characterizing the bulk crystal. 

The Bloch Hamiltonian H{Aa,kz) describes a parity 
and time reversal invariant one dimensional system. In 
the following we will suppress the dependence on A^ and 
consider a purely one dimensional system. To determine 
the end charge N we introduce the "cutting procedure" 
depicted in Fig. [Tfa). We begin with a large but finite 
system with periodic boundary conditions. We then re- 
place the hopping amplitudes ti for all bonds that cross 
the cleavage plane z = by Ati, where A is real. Pro- 
vided z — corresponds to an inversion plane, the one 
dimensional Hamiltonian retains inversion and time re- 
versal symmetry for all A. The fully cleaved crystal cor- 
responds to A = 0. 

For A = 1 the system is translationally invariant, so 
the excess charge near z = 0is(5(A = l) = 0. Since the 
insulator can have no bulk currents, the only way (5(A) 
can change is if a state localized near z — crosses the 
Fermi energy. Thus Q{0) will be the difference between 
the number of states that cross Ep from above and from 
below for A G [0,1]. Kramers' theorem requires that 
every state is at least twofold degenerate, so the number 
of states crossing Ep will be an even integer. Since the 
charge will be divided evenly between the two sides, A^ = 
Q{0)/2 is an integer, which may be written 



A^ = AN: - AN. 



(Al) 



where A A^± is the number of Kramers ' pairs that cross 
Ep from above or below. 

We now relate the parity of A^ to the bulk parity eigen- 
values. To this end it is useful to consider the evolution 
of the spectrum for A G [—1, 1] and to define 



(a) 



. t^^t , 





FIG. 7: (a) A one dimensional inversion symmetric insulator 
cut at 2: = by replacing hopping amplitudes t across z = 
by Xt. The fully cleaved crystal corresponds to A = 0. (b) 
Energy spectrum as a function of A between —1 and 1. The 
conduction and valence bands exchange a Kramers pair of 
states with opposite parity. (c,d) The bulk energy levels at 
A = ±1. For A = — 1 (d) every state at k has a partner at — fc 
with the same energy and opposite parity. For A = -1-1 (c) 
the states at fc = Fi and k = r2 are not paired. 



as the product of the parities of all of the occupied states, 
where each Kramers pair {ijj2a,'4'2a-i) is included only 
once. This quantity is well defined because ^2q = £,2a-i- 
Our proof consists of two steps. We will first show that 



P(1)P(-1) = (-1)^. 
We will then show that 



(A3) 



^(l)^(-l) = 11 [-6m(ri)C2rn(F2)] = 7:. (A4) 



Here C2m(Fi) are the parity of the Bloch states in the 
Tnth Kramers degenerate band at the TRIM k^ = F^, 
and again each Kramers pair is included only once, n^ 
is the number of occupied Kramers degenerate bands. 
Taken together, (JA3| and (JA4| establish the relationship 
summarized by Eqs. I2.1l and l3.5l between the bulk parity 
eigenvalues and the surface fermion parity. 

Eq. IA3I follows from the symmetry of the end state 
spectrum about A = 0. The Hamiltonian H{—X) differs 
H(X) only by a phase twist of tt across z ~ 0. This 
twist can be spread over the entire circumference L by 
performing the gauge transformation 



|^(_A))=e-/^|^(-A)) 



(A5) 



P(A) = n ^2a 

E2a,W<EF 



(A2) 



for < z < L. When L -^ oo the Hamiltonian for 
|V'(— A)) near z — becomes identical to H{\). Thus ev- 
ery bound state lipiW) satisfies Ei{—X) = Ei{X). Since 
(A5) changes the parity, \ipi{X)) and |i/';(~-^)) have op- 
posite parity. 
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It follows that every Kramers pair that crosses the Ep 
at Xq 6 [1,0] has a partner with opposite parity that 
crosses Ep in the opposite direction at —Xq as shown in 
FiglZUb). Thus between A = 1 and A = — 1 the con- 
duction and valence band exchange two Kramers pairs 
with opposite parity, leading to a change in the rela- 
tive sign between P(l) and P(— 1). We conclude that 
P(1)P(-1) = (_l)AW+-HAAr_^ ^jjjpjj Igj^jg directly to 

(US). 

Eq. IA4I follows from a consideration of the parities 
of the Bloch wavefunctions. Consider first the simplest 
case where there is a single Kramers degenerate occupied 
band, as shown in Fig. [7l^c,d). At A = 1 the single 
particle states are labeled by momentum kz — 2fmT/L 
with m = —M/2 + 1, ...,M/2, where M is the number 
of unit cells. At the two TRIM Ti = 0, Ta = Mtt/L 
the parity eigenvalues are ^(Fia). Every other kz has a 



partner — fcz, and even and odd parity combinations of 
the two can be formed. The M/2 — 1 {kz, ~kz) pairs thus 
each contribute —1 to the product in (jA2p . Therefore, 



(A6) 



p(i) = (-i)*^/2-ie(ri)e(r2). 



For A = — 1 the gauge transformation (|A5|) leads to a pe- 
riodic Hamiltonian identical to i?(l), but with momenta 
shifted by 7r/i, as shown in Fig. [7jd). Thus all the 
momenta are paired, so that 



P(-l) = (-l) 



M/2 



(A7) 



Combining (jA6[) and (|A7p leads directly to (jA4|) , which is 
straightforwardly generalized to the case of n;, Kramers 
degenerate bands. 
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